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formed by the Berlin Humboldt university group which has used the exponential relation between 
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1. NSPT and Langevin equation 

It is known that standard diagrammatic Lattice Perturbation Theory (LPT) becomes very com- 
plicated when studying higher orders of typical physical quantities as renormalization factors. 

As an alternative, Numerical Stochastic Perturbation Theory (NSPT) (see e.g. [[!]]) is a powerful 
tool to study higher-loop contributions in LPT: thanks to it, higher-loop results are in fact obtained 
without computing vast numbers of Feynman diagrams. Several applications of NSPT have been 
reported over the last years, for some latest developments see the additional contributions of F. Di 
Renzo, M. Brambilla, C. Torrero and H. Perlt to this conference. 

Here we extend our results reported earlier [§, |||] and study the three-loop ghost propaga- 
tor in Landau gauge to make predictions for standard diagrammatic LPT and compare with non- 
perturbative calculations. 

We use the lattice Langevin equation with stochastic time t 

^-U x ^(t;r 1 )=i(V xs S G [U] + r] x>li {t)) U Xrtl {f,n) , (1.1) 

where r\ is Gaussian random noise, So the gauge action and V x „ the left Lie derivative within 
the gauge group. Discretizing the time t = nx, the equation is integrated numerically in the Euler 
scheme by iteration: 

£/*, M (n+l;T]) = exp(-F JCjjU [£/,?]]) U x ^(n;r]) (1.2) 

with the force 

F X ^[U,T}] =i(TV x ^S G [U} + V^Vx,n)- (1-3) 
Rescaling e = j8t and expanding the gauge links (g oc j3~ 1//2 ) 

/>0 

the Langevin equation at finite time step e turns into a system of updates for each perturbative order 
U x l . The algebra- valued gauge potentials A x41 are related to the gauge lattice link fields U x ijU by 

A x ^=logU x , fl . (1.5) 

Their expansion is given in the form 

Av^Lr'^ (i.6) 

/>o 

and is used to enforce unitarity to all orders in 1 / . 

Each simultaneous Langevin update is augmented by a stochastic gauge-fixing step and by 
subtracting zero modes from A®. From the resulting fields the Green functions of interest can be 
numerically constructed order by order. 

To measure gauge-dependent quantities, exact gauge fixing is needed. We use the Landau 
gauge which is reached by iterative Fourier- accelerated gauge transformations [Q]. 
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2. The ghost propagator in perturbation theory 

2.1 The ghost propagator in NSPT 

It is known that the ghost propagator is defined from the inverse of the Faddeev-Popov (FP) 
operator M. In Landau gauge this operator is constructed by using the lattice covariant D(U) and 
left partial derivatives, M = — V • D(U). Following [||] we use here a definition of M which is most 
suitable for NSPT. 

We introduce the physical lattice momenta p^k^) = ^sin^^-^ = ^sin(^-) and define 
the color diagonal propagator in momentum space as the color trace in the adjoint representation 

(N c = 3) 

In ( 2. 1 ) M~ 1 (k) is the Fourier transform of the inverse FP operator in lattice coordinate space. 



A perturbative expansion is based on the mapping — > — > which allows to 

calculate the inverse FP operator in NSPT recursively (i.e., to any finite order inverting the matrix 
M results in a closed form) 



[M- 1 ] t0) = [mW] ' , [M- 1 ] W = - [m(°)1 ' £ M^-V [M- 1 ] U) . (2.2) 



7=0 

The momentum space ghost propagator is obtained by sandwiching ^ between plane- 

wave vectors. The propagator has to be computed from scratch for each chosen momentum tuple 
(k\ , &2 j kj, , fcj) and different colors of the plane wave. Since we are interested in finding the momen- 
tum behavior as good as possible, these measurements become relatively expensive. 

Multiplying the measured lattice momentum ghost propagator either with (ap) 2 or p 2 , two 
forms of the so-called ghost dressing function are defined: 

7« = M 2 G«, /W=/3 2 G«. (2.3) 

The perturbative construction of M in terms of A differs from the definition adopted in most 
Monte Carlo calculations where a linear relation between the gauge links and gauge potentials is 
used. 

2.2 The ghost propagator in standard LPT 

In the RF-MOM scheme, the renormalized ghost dressing function / RI is defined as 

^W,« M o= /^;,°y\ (2.4) 

with the renormalization condition J Rl ' (p,n, a^y)\ p 2 2 = 1. Therefore, the ghost dressing func- 
tion J(a,p,Ofai) is just the ghost wave function renormalization constant Z(a, J u,a R i/) at A* 2 = q 2 - 
We represent the expansion of the bare J(a,p, CUrf) to n-loop order by 

/ n - lo °P(a,/>,a RI = l + £<4, £ (\&) , if = logM 2 . (2.5) 

i=i k=o V 2 / 
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The leading log coefficients zff coincide with continuum perturbation theory (PT), the subleading 
log coefficients zf{\i>k>o have to be determined both from continuum PT and LPT in the used 
scheme, zfl has to be found in LPT Restricting ourselves in these proceedings to two-loop order 
in Landau gauge and to quenched approximation (compare e.g. we have 

jtf_ 9 m ,_ 315 RI ,_ 2439 35 m , 

The coefficient z^q = 13.8257 has been calculated in [J^J, a first rough estimate of zfg has been 
presented by us at Lattice08 [g]. 

From a RI / = a + (— 221og(aju) +73.9355) 0$ + . .. with the bare coupling Ofo = 

3/(87T 2 j8), we get for the dressing function 

J 2 - l0 °P(a,p,p) = l + i(7 M Jf+7 1)0 ) + ^(7 2i2 j5f 2 +/ 2il j5f+/ 2i0 ) (2.7) 

with 

=-0.0854897, / 1(0 = 0.525314, J 2 ,% = 0.0215195 , J 2 ,i = -0.358423 . (2.8) 

Among others we have to confirm /^o from standard LPT and want to find a precise number for 
/ 2 in Landau gauge quenched QCD using NSPT. 



3. Results 



3.1 Practice of measurements 

Very precise measurements for different lattice sizes and different Langevin steps e are needed. 
Typically we have of the order of 1000 measurements for each momentum tuple. Already at finite e 
the non-integer n = 1/2 (no-loop) contributions to the dressing function have to become negligible 
compared to the neighboring loop contributions. Examples of the dressing function / for n = 1,2,3 
and n = 3/2 vs. p 2 at different volumes and e = 0.01 are shown in Figure [lj 

20 
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f 

Figure 1: Measured ghost dressing function J vs. p 2 for all inequivalent lattice momentum 4-tuples near 
diagonalforL = 8,10,12,14,16,20 and e= 0.01. Left: The one-loop (/3" 1 ), two-loop (/3~ 2 ) and three-loop 
(j3~ 3 ) conttibutions, right: the vanishing (« /3~ 3//2 ) conttibution. 
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We have to take the zero Langevin step limit e — > for each 4-tuple (^1,^2^3,^4) from the 
available finite e measurements at fixed lattice size. This is shown for a particular momentum tuple 
at L = 16 in Figure ||. In order to make contact with standard LPT, the limits L — > 00 and ap — ► 
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Figure 2: Extrapolation to e = of one- (left) and two-loop (right) ghost dressing function for lattice size 
1 6 4 and momentum tuple (1,1,1,1): the fitting function contains both a linear and a quadratic term in £ . 

have to be performed additionally. 

3.2 Fitting logarithmic quantities on finite L 

The dressing function J at zero Langevin step still suffers from finite 0(ap) and 0(pL) correc- 
tions. Our aim is to extract the finite constants in the power expansion of the lattice ghost dressing 
function /, o with very high accuracy. In [ |10| ] it was pointed out that finite-size effects can be 
large when an anomalous dimension comes into play. Having at hand a variety of lattice sizes, we 
address a careful assessment of these effects (the main ideas entering the procedure can be found 



in Ql lfl). Without entering into details, we summarize that strategy of fitting simultaneously 0(ap) 
and O(pL) corrections together using several lattice sizes. 

First we subtract all logarithmic pieces (supposed to be universal and known) from the dressing 
function for each momentum tuple and all lattice sizes. Next we select a range in (ap) 2 = L( fl, / 7 /j) 2 > 
Pn = k^{2n/L) with p\^ m < p 2 < Pmw Within that range we identify a set S of momentum tuples 
(k\,k2,k2,k4) which is common to all chosen lattice sizes. The data in that set are assumed to have 
the same pL effects. Since finite-volume effects decrease with increasing momentum squared, we 
choose as reference fitting point - for an assumed behavior at L = oo - the data point at p 2 « p^ ax 
from the largest lattice size at our disposal. 

Next we perform a non-linear fit using all data points of different L from that set S and the ref- 
erence point collecting for finite size (C m , no functional form) and assuming a functional behavior 
for H (4) (p'^^pp: 

J ifl (ki,k 2 ,k 3 ,k 4 ,p 2 ) =J iQ + a 2 

+ £ C m 8[m,{ki,k2,k 3 ,k 4 }] , 

mes 



i 



~ 2 ~P 

up +r— 
p 

1,2,... 



Pipr + ftp^ + d^)^ 



(3.1) 



Finally we vary the momentum squared window and find an optimal % 2 region which allows 
us to find the 'best' In Figure [3] we show such a # 2 -behavior for a non-linear fit to the one-loop 
dressing function. 
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Figure 3: ;£ 2 -behavior of the non-linear fit for the one-loop dressing function in possible lattice momentum 
windows. 
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Figure 4: L = 8, . . . ,20, S = 7; black filled circles: log-subtracted data from allowed tuples; blue stars: pL 
effects removed; red open circles: C, = 0, H(A) effects removed. Left: one-loop order; right: two-loop order. 

An example of a combined fit at low # 2 is shown in Figure |] for the first and the second loop 
as a function of (ap) 2 . We observe that the numerical data at s = from the chosen set - with all 
logarithmic contributions subtracted - scatter significantly (black filled circles). Switching off the 
O(pL) corrections (C m = in the fit form (fO|)), the blue stars line up in 'rows' according to the 
different hypercubic invariants at infinite lattice volume. The reference point (here the rightmost 
point) is of course unchanged. Finally, after removing also the non-rotational hypercubic effects 
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(leaving only a ^ in (3.1)), we obtain a smooth (almost linear) curve formed by the red open 
circles which directly points towards the fitting constant Jj Q in the zero lattice spacing limit. 

To make a more realistic estimate of 7, o, we have collected fit results from five different sets S 
with minimal % 2 . From these sets we obtain the following (preliminary) constants in Landau gauge 
(Jffi* = 0.525314) 

7 1|0 = 0.52520(46) , 7 2 ,o = 1.489(4) . (3.2) 

The two-loop constant can be transformed into the RF-MOM scheme in that gauge. Our prediction 
for that small contribution is dpo = 9.2(2.7). 

Using the found J2.0 as input for the non-leading log contribution to three loops, we are in the 
position to estimate 73 as well. This analysis is in progress. 
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3.3 Comparison to Monte Carlo data 



Using the A = — log?/ definition as in NSPT, the Berlin Humboldt university group has pro- 
duced Monte Carlo results for the ghost and gluon propagator in Landau gauge and different gauge 



couplings [12]. Since it is assumed that non-perturbative contributions dominate mainly the in- 
frared, it is of interest to compare directly the perturbative ghost dressing function obtained in 
NSPT with its Monte Carlo counterpart for each common momentum tuple. 

We calculate the perturbative dressing function at a given lattice volume summed up to loop 
order w max for a given lattice coupling /3 as follows: 



1 



y — /w 



(3.3) 



In Figure || we compare the perturbative ghost dressing function at lattice size L = 16 to 
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Figure 5: Three-loop NSPT of the ghost dressing function in comparison to Monte Carlo at two j3 values. 

Monte Carlo data at two different /3 values. We observe that in not less than three-loop accuracy 
the perturbative ghost propagator at larger /3 is approximately able to describe the full two-point 
function in the large momentum squared region. 

The situation becomes even worse in comparison with Monte Carlo when trying to define a 
perturbative running coupling from 



a: 



3-loop 



(3-4) 



Here G is the gluon dressing function (see Jl3| , used in NSPT in the same accuracy). This is 
demonstrated in Figure §. So more loops would be necessary to find a satisfactory agreement with 
the non-perturbative data at largest lattice momenta. 



4. Summary 

We have presented a detailed perturbative calculation of the lattice ghost propagator in Landau 
gauge using NSPT. The one-loop constant J\ q perfectly agrees with known V — > °° result. The two- 
loop constant 72,0 is determined with good accuracy for the first time. We have performed a very 



7 



The lattice ghost propagator in Landau gauge up to three loops using NSPT 



A. Schiller 



CO 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 





i 3- loop NSPT 




1 Berlin MO. selected P i— . 




' Berlin MC • 




12 4 , (3 = 6 




i ^ 


□ 


* ^ 



0.12 






3-loopNSPT «^ 
Berlin MC. selected P . — ■ — . - 
' Berlin MC ■ - 


0.11 






12\ (3 = 9 


0.1 






i I 


0.09 


t 


1 


t • 


0.08 






■ » * I Is 


0.07 









10 



P 2 



1(1 



Figure 6: Three-loop NSPT of a s in comparison to Monte Carlo (partly with specially chosen Polyakov 
sectors) for two j3 values. 

careful analysis of all necessary limits. A technique to simultaneously deal with both 0{ap) and 
0{pL) corrections is described in some detail. A comparison with Monte Carlo data of the ghost 
propagator and the running coupling shows that additional loops are needed to better describe the 
asymptotically prevailing perturbative tail at large lattice momenta. 
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